***********************************************************************
** Arctic_code.do ********* January 16, 2025 **************************
** 
** Olczak, Nicholas, & Lisa Dellmuth. 2025. 
** Elite concerns and attitudes about problem solving through
** international organizations. Global Studies Quarterly, 2025.
**
** Input: 
** 		Arctic_survey.dta
**
** Machine: Macbook Pro, Sonoma 14.5
** Program Version: Stata 18.5
**
************************************************************************
************************************************************************ 

ssc install asdoc
net install http://www.stata-journal.com/software/sj18-3/gr0073/
ssc install grstyle
ssc install palettes

* open data

use Arctic_survey.dta, clear

* recode

* importance
tab q0010_0001 
recode q0010_0001 5=. 4=3 3=2 2=1 1=0
recode q0010_0002 5=. 4=3 3=2 2=1 1=0
recode q0010_0003 5=. 4=3 3=2 2=1 1=0
recode q0010_0004 5=. 4=3 3=2 2=1 1=0
recode q0010_0005 5=. 4=3 3=2 2=1 1=0
recode q0010_0008 5=. 4=3 3=2 2=1 1=0 

* confidence
tab q0014_0001
tab q0014_0001, nol
recode q0014_0001 5=. 4=3 3=2 2=1 1=0
gen confgov = q0014_0001
recode q0014_0002 5=. 4=3 3=2 2=1 1=0
gen confEU = q0014_0002
recode q0014_0003 5=. 4=3 3=2 2=1 1=0
gen confUN = q0014_0003
recode q0014_0004 5=. 4=3 3=2 2=1 1=0
gen confAC = q0014_0004
recode q0014_0005 5=. 4=3 3=2 2=1 1=0 
gen confNATO= q0014_0005

gen impdAC=.
replace impdAC=0 if q0010_0004==1|q0010_0004==2
replace impdAC=1 if q0010_0004==3|q0010_0004==4

gen impdNATO=.
replace impdNATO=0 if q0010_0005==1|q0010_0005==2
replace impdNATO=1 if q0010_0005==3|q0010_0005==4

gen impdEU=.
replace impdEU=0 if q0010_0002==1|q0010_0002==2
replace impdEU=1 if q0010_0002==3|q0010_0002==4

gen impdUN=.
replace impdUN=0 if q0010_0003==1|q0010_0003==2
replace impdUN=1 if q0010_0003==3|q0010_0003==4

gen impdgov=.
replace impdgov=0 if q0010_0001==1|q0010_0001==2
replace impdgov=1 if q0010_0001==3|q0010_0001==4

gen confACdum=.
replace confACdum=0 if confAC==0|confAC==1
replace confACdum=1 if confAC==2|confAC==3

gen confNATOdum=.
replace confNATOdum=0 if confNATO==0|confNATO==1
replace confNATOdum=1 if confNATO==2|confNATO==3

gen confEUdum=.
replace confEUdum=0 if confEU==0|confEU==1
replace confEUdum=1 if confEU==2|confEU==3

gen confUNdum=.
replace confUNdum=0 if confUN==0|confUN==1
replace confUNdum=1 if confUN==2|confUN==3

gen confgovdum=.
replace confgovdum=0 if confgov==0|confgov==1
replace confgovdum=1 if confgov==2|confgov==3

* concerns
tab q0021_0001
tab q0021_0001, nol
recode q0021_0001 10=9 9=8 8=7 7=6 6=5 5=4 4=3 3=2 2=1 1=0
recode q0021_0002 10=9 9=8 8=7 7=6 6=5 5=4 4=3 3=2 2=1 1=0
recode q0021_0003 10=9 9=8 8=7 7=6 6=5 5=4 4=3 3=2 2=1 1=0
recode q0021_0004 10=9 9=8 8=7 7=6 6=5 5=4 4=3 3=2 2=1 1=0
recode q0021_0005 10=9 9=8 8=7 7=6 6=5 5=4 4=3 3=2 2=1 1=0
recode q0021_0006 10=9 9=8 8=7 7=6 6=5 5=4 4=3 3=2 2=1 1=0
 
rename q0013_0010 nationalinterest
rename q0027_0004 closetoArctic



* correlations (Tables 1 and 2)

* importance gov't

asdoc pwcorr q0010_0001 q0021_0001 , obs sig star(all) replace save(Table1importance.doc) // climate
asdoc pwcorr q0010_0001 q0021_0002 , obs sig star(all) replace save(Table2importance.doc) // econ growth
asdoc pwcorr q0010_0001 q0021_0003, obs sig star(all) replace save(Table3importance.doc) // military conflicts inv Sweden
asdoc pwcorr q0010_0001  q0021_0004, obs sig star(all) replace save(Table4importance.doc) // militarization of Arctic
asdoc pwcorr q0010_0001  q0021_0005, obs sig star(all) replace save(Table5importance.doc) // env
asdoc pwcorr q0010_0001  q0021_0006, obs sig star(all) replace save(Table6importance.doc) // HR

*EU
asdoc pwcorr q0010_0002 q0021_0001 , obs sig star(all) replace save(Table1aimportance.doc)
asdoc pwcorr q0010_0002 q0021_0002 , obs sig star(all) replace save(Table2aimportance.doc)
asdoc pwcorr q0010_0002 q0021_0003, obs sig star(all) replace save(Table3aimportance.doc)
asdoc pwcorr q0010_0002  q0021_0004, obs sig star(all) replace save(Table4aimportance.doc) 
asdoc pwcorr q0010_0002  q0021_0005, obs sig star(all) replace save(Table5aimportance.doc) 
asdoc pwcorr q0010_0002  q0021_0006, obs sig star(all) replace save(Table6aimportance.doc)

*UN
asdoc pwcorr q0010_0003 q0021_0001 , obs sig star(all) replace save(Table1bimportance.doc)
asdoc pwcorr q0010_0003 q0021_0002 , obs sig star(all) replace save(Table2bimportance.doc)
asdoc pwcorr q0010_0003 q0021_0003, obs sig star(all) replace save(Table3bimportance.doc)
asdoc pwcorr q0010_0003  q0021_0004, obs sig star(all) replace save(Table4bimportance.doc) 
asdoc pwcorr q0010_0003  q0021_0005, obs sig star(all) replace save(Table5bimportance.doc) 
asdoc pwcorr q0010_0003  q0021_0006, obs sig star(all) replace save(Table6bimportance.doc)

*AC
asdoc pwcorr q0010_0004 q0021_0001 , obs sig star(all) replace save(Table1cimportance.doc)
asdoc pwcorr q0010_0004 q0021_0002 , obs sig star(all) replace save(Table2cimportance.doc)
asdoc pwcorr q0010_0004 q0021_0003, obs sig star(all) replace save(Table3cimportance.doc)
asdoc pwcorr q0010_0004  q0021_0004, obs sig star(all) replace save(Table4cimportance.doc) 
asdoc pwcorr q0010_0004  q0021_0005, obs sig star(all) replace save(Table5cimportance.doc) 
asdoc pwcorr q0010_0004  q0021_0006, obs sig star(all) replace save(Table6cimportance.doc)

*NATO
asdoc pwcorr q0010_0005 q0021_0001 , obs sig star(all) replace save(Table1dimportance.doc)
asdoc pwcorr q0010_0005 q0021_0002 , obs sig star(all) replace save(Table2dimportance.doc)
asdoc pwcorr q0010_0005 q0021_0003, obs sig star(all) replace save(Table3dimportance.doc)
asdoc pwcorr q0010_0005  q0021_0004, obs sig star(all) replace save(Table4dimportance.doc) 
asdoc pwcorr q0010_0005  q0021_0005, obs sig star(all) replace save(Table5dimportance.doc) 
asdoc pwcorr q0010_0005  q0021_0006, obs sig star(all) replace save(Table6dimportance.doc)

* confidence

*confidence gov't
asdoc pwcorr q0014_0001 q0021_0001 , obs sig star(all) replace save(Table1confidence.doc)
asdoc pwcorr q0014_0001 q0021_0002 , obs sig star(all) replace save(Table2confidence.doc)
asdoc pwcorr q0014_0001 q0021_0003, obs sig star(all) replace save(Table3confidence.doc)
asdoc pwcorr q0014_0001  q0021_0004, obs sig star(all) replace save(Table4confidence.doc) 
asdoc pwcorr q0014_0001  q0021_0005, obs sig star(all) replace save(Table5confidence.doc) 
asdoc pwcorr q0014_0001  q0021_0006, obs sig star(all) replace save(Table6confidence.doc)

*EU
asdoc pwcorr q0014_0002 q0021_0001 , obs sig star(all) replace save(Table1aconfidence.doc)
asdoc pwcorr q0014_0002 q0021_0002 , obs sig star(all) replace save(Table2aconfidence.doc)
asdoc pwcorr q0014_0002 q0021_0003, obs sig star(all) replace save(Table3aconfidence.doc)
asdoc pwcorr q0014_0002  q0021_0004, obs sig star(all) replace save(Table4aconfidence.doc) 
asdoc pwcorr q0014_0002  q0021_0005, obs sig star(all) replace save(Table5aconfidence.doc) 
asdoc pwcorr q0014_0002  q0021_0006, obs sig star(all) replace save(Table6aconfidence.doc)

*UN
asdoc pwcorr q0014_0003 q0021_0001 , obs sig star(all) replace save(Table1bconfidence.doc)
asdoc pwcorr q0014_0003 q0021_0002 , obs sig star(all) replace save(Table2bconfidence.doc)
asdoc pwcorr q0014_0003 q0021_0003, obs sig star(all) replace save(Table3bconfidence.doc)
asdoc pwcorr q0014_0003  q0021_0004, obs sig star(all) replace save(Table4bconfidence.doc) 
asdoc pwcorr q0014_0003  q0021_0005, obs sig star(all) replace save(Table5bconfidence.doc) 
asdoc pwcorr q0014_0003  q0021_0006, obs sig star(all) replace save(Table6bconfidence.doc)

*AC
asdoc pwcorr q0014_0004 q0021_0001 , obs sig star(all) replace save(Table1cconfidence.doc)
asdoc pwcorr q0014_0004 q0021_0002 , obs sig star(all) replace save(Table2cconfidence.doc)
asdoc pwcorr q0014_0004 q0021_0003, obs sig star(all) replace save(Table3cconfidence.doc)
asdoc pwcorr q0014_0004  q0021_0004, obs sig star(all) replace save(Table4cconfidence.doc) 
asdoc pwcorr q0014_0004  q0021_0005, obs sig star(all) replace save(Table5cconfidence.doc) 
asdoc pwcorr q0014_0004  q0021_0006, obs sig star(all) replace save(Table6cconfidence.doc)

*NATO
asdoc pwcorr q0014_0005 q0021_0001 , obs sig star(all) replace save(Table1dconfidence.doc)
asdoc pwcorr q0014_0005 q0021_0002 , obs sig star(all) replace save(Table2dconfidence.doc)
asdoc pwcorr q0014_0005 q0021_0003, obs sig star(all) replace save(Table3dconfidence.doc)
asdoc pwcorr q0014_0005  q0021_0004, obs sig star(all) replace save(Table4dconfidence.doc) 
asdoc pwcorr q0014_0005  q0021_0005, obs sig star(all) replace save(Table5dconfidence.doc) 
asdoc pwcorr q0014_0005  q0021_0006, obs sig star(all) replace save(Table6dconfidence.doc)



*Figure 1

grstyle clear
set scheme plotplain
grstyle init
grstyle set plain, box
grstyle color background white 
grstyle yesno draw_major_hgrid yes
grstyle yesno draw_major_ygrid yes
grstyle color major_grid gs8
grstyle linepattern major_grid dot
grstyle set legend 2, box inside
grstyle color ci_area gs12%50

#delimit ;
graph bar q0021_0003 q0021_0004 q0021_0001 q0021_0002 q0021_0005 q0021_0006, 
 ylabel(0 "0-Not at all worried" 1 "1" 2 "2" 3 "3" 4 "4" 5 "5" 6 "6" 7 "7" 8 "8" 9 "9-Very worried", 
labsize(small) labcolor(black) angle(horizontal) ///
glcolor(gs14) gmax) bar(3, color(g2))  
bar(1, fcolor(gs14) fintensity(100)) 
bar(2, fcolor(gs8) fintensity(100)) 
bar(3, fcolor(gs2) fintensity(100)) 
blabel(bar, format(%9.1f) size(small) lcolor(none))
;
#delimit cr

sum q0021_0003 q0021_0004 q0021_0001 q0021_0002 q0021_0005 q0021_0006
ttest q0021_0003 = q0021_0004    
ttest q0021_0003 =q0021_0001
ttest q0021_0003 =q0021_0002
ttest q0021_0003 =q0021_0005
ttest q0021_0003 =q0021_0006

ttest q0021_0004 =q0021_0001
ttest q0021_0004 =q0021_0002
ttest q0021_0004 =q0021_0005
ttest q0021_0004 =q0021_0006

ttest q0021_0001 =q0021_0002
ttest q0021_0001 =q0021_0005
ttest q0021_0001 =q0021_0006

ttest q0021_0002 =q0021_0006
ttest q0021_0002 =q0021_0005
ttest q0021_0005 =q0021_0006


*Figure 2

#delimit ;
graph bar q0010_0001 q0010_0002 q0010_0003 q0010_0004 q0010_0005, 
  ylabel(0 "0-Not at all important" 1 "1-Not very important" 2 "2-Quite important" 3 "3-Very important", 
labsize(small) labcolor(black) angle(horizontal) ///
glcolor(gs14) gmax) bar(3, color(g2))  
 ytitle(Confidence)
bar(1, fcolor(gs14) fintensity(100)) 
bar(2, fcolor(gs8) fintensity(100)) 
bar(3, fcolor(gs2) fintensity(100)) 
blabel(bar, format(%9.1f) size(small) lcolor(none))
;
#delimit cr

sum q0010_0001 q0010_0002 q0010_0003 q0010_0004 q0010_0005
ttest q0010_0001 =q0010_0002   
ttest q0010_0001 =q0010_0003
ttest q0010_0001 =q0010_0004
ttest q0010_0001 =q0010_0005

ttest q0010_0002 =q0010_0003
ttest q0010_0002 =q0010_0004
ttest q0010_0002 =q0010_0005

ttest q0010_0003 =q0010_0004
ttest q0010_0003 =q0010_0005
ttest q0010_0004 =q0010_0005

*Figure 3

#delimit ;
graph bar q0014_0001 q0014_0002 q0014_0003 q0014_0004 q0014_0005, 
 ylabel(0 "0.None at all" 1 "1-Not very much" 2 "2-Quite much" 3 "3-Very much", 
labsize(small) labcolor(black) angle(horizontal) ///
glcolor(gs14) gmax) bar(3, color(g2))   
blabel(bar, format(%9.1f) size(small) lcolor(none))
;
#delimit cr

sum q0014_0001 q0014_0002 q0014_0003 q0014_0004 q0014_0005
ttest q0014_0001 =q0014_0002   
ttest q0014_0001 =q0014_0003
ttest q0014_0001 =q0014_0004
ttest q0014_0001 =q0014_0005

ttest q0014_0002 =q0014_0003
ttest q0014_0002 =q0014_0004
ttest q0014_0002 =q0014_0005

ttest q0014_0003 =q0014_0004
ttest q0014_0003 =q0014_0005
ttest q0014_0004 =q0014_0005



* ROBUSTNESS CHECKS


* Table E1 

* add: nationalist values, geographic identity, confgovdum, age, gender, edu

* importance

ologit q0010_0004 q0021_0003 q0021_0004 nationalinterest    q0027_0006 confgov gender agecat, robust 
eststo AC  

ologit q0010_0005 q0021_0003 q0021_0004 nationalinterest     q0027_0006  confgov gender agecat, robust
eststo NATO 

ologit q0010_0002 q0021_0003 q0021_0004 nationalinterest     q0027_0006 confgov gender agecat, robust
eststo EU 

ologit q0010_0003 q0021_0003 q0021_0004 nationalinterest    q0027_0006 confgov gender agecat, robust
eststo UN 

esttab AC NATO EU UN  


* confidence 

ologit confAC q0021_0003 q0021_0004 nationalinterest    q0027_0006 confgov gender agecat, robust
eststo ACc   

ologit confNATO q0021_0003 q0021_0004 nationalinterest    q0027_0006 confgov gender agecat, robust
eststo NATOc 

ologit confEU q0021_0003 q0021_0004 nationalinterest     q0027_0006 confgov gender agecat, robust
eststo EUc 

ologit confUN q0021_0003 q0021_0004 nationalinterest     q0027_0006 confgov gender agecat, robust
eststo UNc 
 
esttab ACc NATOc EUc UNc  


#delimit ;
esttab AC NATO EU UN ACc NATOc EUc UNc using ologit.rtf, b(%12.3f)
se(%12.3f) starlevels(* .10 ** .05 *** .01) scalars(ll) bic ar2(%12.3f) 
replace;
#delimit cr



* Table E2 

* importance

ologit q0010_0004 q0021_0003 q0021_0004 q0021_0001 q0021_0002 q0021_0005 q0021_0006 nationalinterest    q0027_0006 q0014_0001 gender agecat, robust 
eststo AC  

ologit q0010_0005 q0021_0003 q0021_0004 q0021_0001 q0021_0002 q0021_0005 q0021_0006 nationalinterest     q0027_0006  q0014_0001 gender agecat, robust
eststo NATO 

ologit q0010_0002 q0021_0003 q0021_0004 q0021_0001 q0021_0002 q0021_0005 q0021_0006 nationalinterest     q0027_0006 q0014_0001 gender agecat, robust
eststo EU 

ologit q0010_0003 q0021_0003 q0021_0004 q0021_0001 q0021_0002 q0021_0005 q0021_0006 nationalinterest    q0027_0006 q0014_0001 gender agecat, robust
eststo UN 

esttab AC NATO EU UN  

* confidence

ologit confAC q0021_0003 q0021_0004 q0021_0001 q0021_0002 q0021_0005 q0021_0006 nationalinterest    q0027_0006 confgov gender agecat, robust
eststo ACc   

ologit confNATO q0021_0003 q0021_0004 q0021_0001 q0021_0002 q0021_0005 q0021_0006 nationalinterest    q0027_0006 confgov gender agecat, robust
eststo NATOc 

ologit confEU q0021_0003 q0021_0004 q0021_0001 q0021_0002 q0021_0005 q0021_0006 nationalinterest     q0027_0006 confgov gender agecat, robust
eststo EUc 

ologit confUN q0021_0003 q0021_0004 q0021_0001 q0021_0002 q0021_0005 q0021_0006 nationalinterest     q0027_0006 confgov gender agecat, robust
eststo UNc 
 
esttab ACc NATOc EUc UNc  


#delimit ;
esttab AC NATO EU UN ACc NATOc EUc UNc using ologit_wadditionalconcerns.rtf, b(%12.3f)
se(%12.3f) starlevels(* .10 ** .05 *** .01) scalars(ll) bic ar2(%12.3f) 
replace;
#delimit cr


* Table E3

* add: elite categories
tab cat
tab cat, nol
gen cat4=.
replace cat4=1 if cat==1 // CS
replace cat4=2 if cat==2 // bus
replace cat4=3 if cat==3|cat==4 // pol and admin
replace cat4=4 if cat==4 // research

* importance

ologit q0010_0004 q0021_0003  q0021_0004 i.cat4 nationalinterest    q0027_0006 confgov gender agecat, robust 
eststo AC  

ologit q0010_0005 q0021_0003   q0021_0004 i.cat4 nationalinterest     q0027_0006  confgov gender agecat, robust
eststo NATO 

ologit q0010_0002 q0021_0003  q0021_0004 i.cat4  nationalinterest     q0027_0006 confgov gender agecat, robust
eststo EU 

ologit q0010_0003 q0021_0003  q0021_0004 i.cat4 nationalinterest    q0027_0006 confgov gender agecat, robust
eststo UN 

esttab AC NATO EU UN  

* confidence

ologit confAC q0021_0003 q0021_0004 i.cat4 nationalinterest    q0027_0006 confgov gender agecat, robust
eststo ACc   

ologit confNATO q0021_0003 q0021_0004 i.cat4 nationalinterest    q0027_0006 confgov gender agecat, robust
eststo NATOc 

ologit confEU q0021_0003 q0021_0004 i.cat4 nationalinterest     q0027_0006 confgov gender agecat, robust
eststo EUc 

ologit confUN q0021_0003 q0021_0004 i.cat4 nationalinterest     q0027_0006 confgov gender agecat, robust
eststo UNc 
 
esttab ACc NATOc EUc UNc  


#delimit ;
esttab AC NATO EU UN ACc NATOc EUc UNc using ologit_cat4.rtf, b(%12.3f)
se(%12.3f) starlevels(* .10 ** .05 *** .01) scalars(ll) bic ar2(%12.3f) 
replace;
#delimit cr


* Table E4 

* add: elite cat with interaction term with national security worries

gen cat2=.
replace cat2= 0 if cat==1|cat==2|cat==4 // societal elites
replace cat2= 1 if cat==3|cat==4 // political elites


* importance

ologit q0010_0004 c.q0021_0003##i.cat2 q0021_0004 nationalinterest q0027_0006 confgov gender agecat, robust
eststo AC  

ologit q0010_0005 c.q0021_0003##i.cat2 q0021_0004 nationalinterest     q0027_0006  confgov gender agecat, robust
eststo NATO 

ologit q0010_0002 c.q0021_0003##i.cat2 q0021_0004 nationalinterest     q0027_0006 confgov gender agecat, robust
eststo EU 

ologit q0010_0003 c.q0021_0003##i.cat2 q0021_0004 nationalinterest    q0027_0006 confgov gender agecat, robust
eststo UN 

esttab AC NATO EU UN  

* confidence

ologit confAC c.q0021_0003##i.cat2 q0021_0004 nationalinterest    q0027_0006 confgov gender agecat, robust
eststo ACc   

ologit confNATO c.q0021_0003##i.cat2 q0021_0004 nationalinterest    q0027_0006 confgov gender agecat, robust
eststo NATOc 

ologit confEU c.q0021_0003##i.cat2 q0021_0004 nationalinterest     q0027_0006 confgov gender agecat, robust
eststo EUc 

ologit confUN c.q0021_0003##i.cat2 q0021_0004 nationalinterest     q0027_0006 confgov gender agecat, robust
eststo UNc 
 
esttab ACc NATOc EUc UNc  


#delimit ;
esttab AC NATO EU UN ACc NATOc EUc UNc using ologit_int.rtf, b(%12.3f)
se(%12.3f) starlevels(* .10 ** .05 *** .01) scalars(ll) bic ar2(%12.3f) 
replace;
#delimit cr




* Table E5

* add: elite cat with interaction term with general security worries


ologit q0010_0004 c.q0021_0004##i.cat2 q0021_0003 nationalinterest    q0027_0006 confgov gender agecat, robust 
eststo AC  

ologit q0010_0005 c.q0021_0004##i.cat2 q0021_0003 nationalinterest     q0027_0006  confgov gender agecat, robust
eststo NATO 

ologit q0010_0002 c.q0021_0004##i.cat2 q0021_0003 nationalinterest     q0027_0006 confgov gender agecat, robust
eststo EU 

ologit q0010_0003 c.q0021_0004##i.cat2 q0021_0003 nationalinterest    q0027_0006 confgov gender agecat, robust
eststo UN 

esttab AC NATO EU UN  

* confidence

ologit confAC c.q0021_0004##i.cat2 q0021_0003 nationalinterest    q0027_0006 confgov gender agecat, robust
eststo ACc   

ologit confNATO c.q0021_0004##i.cat2 q0021_0003 nationalinterest    q0027_0006 confgov gender agecat, robust
eststo NATOc 

ologit confEU c.q0021_0004##i.cat2 q0021_0003 nationalinterest     q0027_0006 confgov gender agecat, robust
eststo EUc 

ologit confUN c.q0021_0004##i.cat2 q0021_0003 nationalinterest     q0027_0006 confgov gender agecat, robust
eststo UNc 
 
esttab ACc NATOc EUc UNc  


#delimit ;
esttab AC NATO EU UN ACc NATOc EUc UNc using ologit_intwgensec.rtf, b(%12.3f)
se(%12.3f) starlevels(* .10 ** .05 *** .01) scalars(ll) bic ar2(%12.3f) 
replace;
#delimit cr




